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ABSTRACT 

In certain two-dimensional models, collapsing matter forms a black hole if and 
only if the incoming energy flux exceeds the Hawking radiation rate. Near the 
critical threshold, the black hole mass is given by a universal formula in terms 
of the distance from criticality, and there exists a scaling solution describing the 
formation and evaporation of an arbitrarily small black hole. 
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Recently remarkable critical behavior has been discovered at the onset of clas- 
sical black hole formation. Choptuik [1] considered a two-dimensional theory ob- 
tained from the S'-wave sector of four- dimensional general relativity coupled to a 
massless scalar field. A sufficiently weak incoming S'-wave pulse is simply reflected 
through the origin to an outgoing pulse. This behavior changes qualitatively as 
the initial amplitude is increased. Above a certain threshold, the pulse crosses its 
own Schwarzschild radius before it reaches the origin, and a black hole is formed. 
Choptuik [1] numerically investigated the mass Mbh of the resulting black hole as 
a function of the distance S in the initial data space from the threshold. For small 
S he finds: 

logMBH = llog6 + 0{6^), (1) 

where the critical exponent 7 is numerically found to be near .37. This result 
appears quite universal and is insensitive to the precise definition of S or scalar 
field couplings. In addition the near-threshold scaling solution was found to have 
fascinating self-similar oscillations. Related work can be found in [2]. At present 
there appears to be little analytic or conceptual understanding of these interesting 
phenomena. 

Two-dimensional theories obtained by S-wave reduction from four dimensions 
have also been of recent interest as simplified arenas for the study of quantum black 
hole evaporation [3]. In this context null matter which reduces to a free conformal 
field in two dimensions is usually considered because it is much simpler than a 
scalar field (which reduces to a field with complicated gravitational couplings) . At 
the classical level, every incoming S-wave pulse of such null matter, no matter how 
weak, forms a black hole. A threshold appears, however, when quantum effects 
are incorporated: Energy must be thrown at the origin at a rate faster than a 
newly-formed black hole wants to Hawking evaporate. 

In this paper we study the onset of quantum black hole formation using the 
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semiclassically soluble two-dimensional RST model [4]. We find analytically that 

log ^ = 1 log 5 - i log a + 0{5^) . (2) 

Mbh is defined here as the incoming energy of the null matter swallowed by the 
black hole during its lifetime and A is the dimensionful parameter of the model. The 
0{5^) term depends universally on a, the second derivative of the energy density 
at the point where the critical threshold is exceeded*. The 0{5^/'^) corrections 
depend non-universally on the shape of the incoming pulse. We also find a scaling 
solution near criticality which corresponds to the formation and evaporation of an 
arbitrarily small black hole. We have not determined whether the scaling (2) is 
universal with respect to small changes in the coupling constants of the theory. 
Presumably numerical work is required to answer this interesting question. 

There are obvious similarities between our results and those of Ref. [1], but 
there are also apparent differences. First of all, our critical exponent is a rational 
number whereas Choptuik's at least appears to be irrational. Second, there is 
no analog of the self-similar oscillations in our work. This could be a special 
feature arising from the linear nature of the RST equations, while more general 
two-dimensional models (which are not analytically soluble) might exhibit such 
oscillations. 

We now present a derivation of the scaling relation (2). The semi-classical 
effective action for the RST model ^ is 



S=— d^x 
27r 



^ i=\ (3) 
~ 9^/ ^^^V^-^ / d^x'V-gi^') R{x)G{x-x')R{x') , 

where g^i, is the two-dimensional metric, is a scalar field called the dilaton, fi 
are A^" minimally coupled scalar matter fields and G is a Green function for the 



-k It would be interesting to determine if the 0{5^) term of (1) has a similar universal 
dependence. 

f For more details on this model consult [4,6]. We use the conventions of [6]. 
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operator V^. The effective action includes the one-loop Liouville term due to the 
matter fields and if is large this term provides the dominant quantum back- 
reaction on the geometry. This model differs from the original CGHS model [7] 
by a finite local counterterm, which restores a global symmetry of the classical 
theory and enables writing down exact semi-classical solutions in a rather simple 
form. Numerical analyses of the original model [8] indicate, however, that the two 
models are similar, i.e. that the qualitative behavior of semi-classical solutions is 
not sensitive to the existence of the global symmetry. 

It is convenient to work in conformal gauge, = —^e^P, g++ — g = 0, 

and use the global symmetry to further fix the coordinates to "Kruskal gauge", 
where p — (f)+ ^ log ^. This ehminates the conformal factor from the discussion. 
If we define a new dilaton field, 

Q, = — e + + - log — 

the semi-classical equations reduce to 

d+d^Q = - , (5) 

-dln=Tl^ + t±, 

where — Y^^=i{^±fi)'^ ■ The function takes the value tj^ — —4(^^)2 in 

Kruskal coordinates for any incoming matter energy fiux which vanishes sufficiently 
rapidly at asymptotic early and late times. The field redefinition (4) is degenerate 
— \ and ^ < \ does not correspond to a real value of 0. The curve ^ — \ 
is the analog of the origin of radial coordinates in higher dimensional gravity''' and 
solutions should not be continued beyond it. Instead, RST impose the following 

I This analogy can be made precise when the model is interpreted as an effective theory for 
radial modes of near extreme magnetic dilaton black holes in four dimensional gravity [3] . 
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boundary conditions at this curve, wherever it is timehke [4], 



= . (6) 
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Incoming energy flux is then reflected off the boundary and the total outgoing 
flux (including the anomalous part t_) can be determined by using (6). The 
boundary curve undergoes dynamical motion in response to the incoming matter 
and this gives rise to some non-trivial behavior in this model, in spite of the extreme 
simplicity of the field equations (5). 

The boundary conditions (6) ensure semi-classical energy conservation and also 
that the physical curvature remains flnite at the boundary curve as long as it is 
timelike. It should be noted, however, that these boundary conditions are not 
the most general ones allowed, and our results may depend qualitatively on this 
choice. Indeed (6) is incompatible with Dirichlet or Neumann boundary conditons 
on the fi, and may not be realizable as the semiclassical limit of any fully quantum 
mechanical boundary conditions [5]. Alternate possibilities are currently being 
explored. 

The solution corresponding to incoming matter energy flux, which tapers off 
at early and late times, but is otherwise quite general, is given by 



Q{x+,x-) = -x+{X^x- + P+{x+)) + \M{x+)-\M{x+{x-))-]\og (^%) , 

A ^ ' A ^ '4 x^{x ) 

(7) 

where 

M{x+) ^xj duuTl^{u) , = j duTl^{u) , (8) 



and x'^{x ) is the x'^ value of the point on the boundary curve from which the 
reflected signal propagates to {x'^,x~). Wherever the boundary curve is timelike 
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its shape is a simple function of the incoming energy flux, 



A^xo = - -— . (9) 
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If, however, the incoming energy flux becomes larger than the outgoing Hawking 
flux of a two-dimensional black hole, T:l^{x~^) > j^^^+yi} for some value of x'^ then 
the boundary curve becomes spacelike and it is less trivial to determine its shape. 

Spacelike segments of the boundary curve are curvature singularities. They 
form inside regions of future trapped points which are bounded by an apparent 
horizon, located where d+Q — [9]. Applying d+ to (7) one flnds that the appar- 
ent horizon curve {xjj,x]j) satisfles (9) for all values of xj. In other words, the 
apparent horizon coincides with the boundary curve where the latter is timelike, 
but where the boundary becomes spacelike the two curves separate and the appar- 
ent horizon cloaks the singularity, as shown in figure 1. Once the incoming energy 
flux falls below the threshold value the black hole evaporates and the apparent 
horizon approaches the singularity. The curves join again at the endpoint of the 
evaporation, which is denoted by E in figure 1. The null line-segment x~ = x^, 
x~^{x^) < x~^ < x~^, is the global event horizon of the geometry. We define the 
black hole mass to be the total energy (as measured at I~) which enters the black 
hole 

Mbh = M{x+) - M{x+{xe)) . (10) 

Since both ends of the global horizon are on the boundary curve and also on the 
apparent horizon it immediately follows from (7) and (9) that 

Now consider black hole formation just above threshold. For concreteness 



The black hole temperature is independent of mass in this model and an energy flux of 
tI+{x+) = l/4(x+)2 in Kruskal coordinates corresponds to a uniform incoming flux, 
T:f_l_((7+) — in the coordinate, Act+ = logAa;+, appropriate to asymptotic inertial 

observers. 
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assume that the incoming energy flux has a maximum at Xa^ = logAx^ = 0, 
where its value is t{_|_ = A^(| + 5), and that the flux is below threshold everywhere 
except near cr"*" = so that only a single small black hole is formed. In the scaling 
limit, 5 ^ 0, a. generic incoming flux of this type may be parametrized as follows 
near a'^ — 0, 

TUa+) = X\\ + S){l-a AV+' + ...)■ (12) 

f 

It turns out that the higher terms in the Taylor expansion of Tjf.^ do not contribute 
in the scaling limit. 




FIGURE 1. Kruskal diagram for black hole formation and evaporation in the scaling 
regime. A near-critical flux of matter energy is incident from = — OO. The solid 
curve is the Q — j boundary and the dashed curve is the apparent horizon. 
The spacelike portion of the solid curve is the black hole singularity. 
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To obtain the shape of the apparent horizon curve, we transform to Kruskal 
coordinates, compute P-^-{x~^) for this flux distribution, and insert the result into 
(9). We are interested in the 6 << 1 limit and for that purpose it is convenient to 
shift and rescale the Kruskal coordinates as follows, 

Ax+ = l + y^(a+-2), Ax- = -ip+(i)-i + A/?(a- + |). (13) 

\ a A A A \ a 6 

The origin of the {a'^,a~) coordinate system has been chosen where, to leading 
order in 5, the boundary curve turns spacelike, as shown in figure 1. In the scaling 
region where higher order terms in 6 can be dropped, the apparent horizon takes 
a simple and universal form, 

a]j{a+) = -^a+' + ^af . (14) 

The only reference made to the parameter a is in the definition of the scaling 
variables (13). Higher orders in the Taylor expansion of T^^{a'^) in (12) only 
contribute terms carrying positive powers of S, which can be ignored in the scaling 
region. 

It is straightforward to show that the singularity curve is also universal in this 
limit. Expressing (7) in terms of (a"'",a~) coordinates and applying to both 
sides yields the following differential equation for the — ^ curve, 

diQj 

(a+ - a+(a^)) ^ = -a^ + ajj{a+) . (15) 

This is a non-linear differential equation and we do not have an analytic solution, 
but there is no explicit dependence on the flux parameters so the shape of the curve 
must be universal. However, we do not need the whole singularity curve in order 
to determine the black hole mass (11). It is sufficient to locate the endpoint of 
evaporation and this can be achieved by the following geometric argument, which 
is illustrated in figure 2. 
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FIGURE 2. The two shaded regions have equal area for any point (a^ , a^) on the singularity curve. 

Take any point on the singularity curve and consider the past-directed null 
line = which extends from a+ = to the timclikc boundary curve at 
a"*" = a^(a^). Denote by P the point where this null line intersects the apparent 
horizon. The value of ri{a^,ap) can be obtained by integrating d+Vl along this 
null line starting either from the spacelike singularity at or from the timelike 
boundary curve at a^(a^). The value of d+fl on a~ — can in turn be obtained 
by integrating d-d+ft along null lines of constant a'^ starting from the apparent 
horizon curve (where d+il vanishes by definition). By using the equation of motion, 
= — A^, inside the double integral one thus finds that fl{ap, ap) — | equals 
the area of each of the shaded regions in figure 2 and therefore these areas must 
be equal. Interestingly, this "equal area rule" is an exact relationship which also 
holds outside of the scaling limit. 

The equal area rule, when applied at the global event horizon, a~ — a^, gives 
(a^) . Both endpoints of the global event horizon are also points 
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on the apparent horizon curve (14) so we have — + J2^e ~ ^I^b^'^e)'^ 
^o^(a^)^. From these two relations it follows that — 2 + \/l2, o^(a^) = 
2 — -\/l2, and the black hole mass (11) is to leading order seen to be 



Mbh 
A 




(16) 



which is the scaling relation (2). 

Finally, it is straightforward to integrate (15) numerically to obtain the entire 
singularity curve. The result is plotted in figure 3. 



a minus 
1 




FIGURE 3. Singularity curve obtained from numerical calculation. 
The dashed curve is the apparent horizon (14). 
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